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For many driven-nonequilibrium systems, the probability distribution functions of magnitude and 
recurrence-time of large events follow a powerlaw indicating a strong temporal correlation. In this 
paper we argue why these probability distribution functions are ubiquitous in driven nonequilibrium 
systems, and we derive universal scaling laws connecting the magnitudes, recurrence-time, and 
spatial intervals of large events. The relationships between the scaling exponents have also been 
studied. We show that the ion-channel current in Voltage-dependent Anion Channels obeys the 
universal scaling law. 
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The physics of nonequilibrium systems has been of major interest in last several decades. Recently large events 
f*^ ' have become a major resource for characterizing nonequilibrium system. Bak et al. [J, |2( constructed a model called 
"sandpile" to abstract some of the essential features of driven-nonequilibrium systems. Original model of Bak et 
al. (BTW) as well as its variants exhibit powerlaw behaviour for the magnitude and duration of avalanches, which 
play the role of large events. Later scientists started investigating the temporal correlations in various systems using 
recurrence-time pdf of large events. In the BTW model, the pdf of recurrence time (r) of large events is exp(— t/t c )/t c , 
which indicates that the triggering of large events in this system are uncorrelated |3J. |j|. This is a generic feature of 
Poisson process. It was however soon realized that the events are correlated in many nonequilibrium systems. Boffatta 
et al. 0| studied solar flares, and found that the pdf of recurrence time for the solar flares is a powerlaw, i. e., 

^; p(t) = At-p, (i) 

O 

with p rs 2.4. These authors also studied the recurrence-time distribution for the dissipation-bursts in magnetohy- 
drodynamic turbulence, and again found a powerlaw with j3 ~ 2.7. These results show that the large events in solar 
flares and magnetohydrodynamic turbulence are temporally correlated. Note that the pdf of solar-flare magnitudes is 
a powerlaw just like sandpile model 0,0. The difference between the sandpile model and solar flare is in the temporal 

£h ' correlations between the large events 0, U| ■ 
'r^ , Temporal correlations described above have been observed for many other nonequilibrium systems. Corral |7j, la] 

studied the pdf of recurrence time for large earthquakes (beyond 2.5 in Richter scale), and found it to be a powerlaw 

J> , with (i rs 0.3. Bak et al. J2( had done similar analysis. Earlier Guttenburg and Richter [lfj showed that the magnitude 
of earthquakes follows a powerlaw. "Extremal models" (see Paczuski et al. [ll| for review of these models) which 
are used to describe interface depinning, biological evolution etc. also have similar temporal correlations. Recently 
Banerjee et al. [l2| studied ion-channel currents in Voltage Dependent Anion Channel (VDAC). Here too, the pdf 
of magnitude as well as that of recurrence-time follow a powerlaw with /3 rs 1.5. In thermal convection of Helium, 
Sreenivasan et al. [13] observed large-scale wind which switches directions; the pdf of inter-switch interval is given by 
(1/t) exp(— t/tm). This feature appears to indicate "release of convective stresses" during switching, just like energy 
release in earthquakes. The pdf of recurrence time is connected to the power spectrum of the signal. Several research 
groups [ll|, LL3, LLJ, tLa] have derived 1// power spectrum using the powerlaw distribution of recurrence-time. 

There is a recent work on spatial correlations in driven nonequilibrium systems. Davidsen and Paczuski [la] 
computed the correlation between the epicenters of subsequent earthquakes, and found the pdf for the distance Ar 
between the subsequent epicenters to be a powerlaw. Earlier, Krishnamurthy and Barma |l7t Il8l | found a powerlaw 
distribution for the long-range jumps in self-organized interface depinning model. 



Corral [7, £j and Bak et al. |2| combined the pdfs of magnitude and recurrence-time and derived a unified scaling 
law. For events greater than a threshold I c , then the pdf for recurrence-time is 

P{r) = Rf(rR) (2) 

with 

/T \ -0.33 

J{tR) a OR- 033 (-) exp(-ri?/a) (3) 

and R — RolO~ bMc (b rs 1) is the average number of earthquakes (with I > I c ) per unit time, and a is a constant. 11 
these work however are model specific. In our paper we present a general scaling laws which should be applicable to 
all the above systems and other driven nonequilibrium systems. Here we also derive universal relationships between 
the exponents. Note that M c — log 10 A c is the Richter scale. 

Eqs. <|2l EJ form the unified scaling law for the earthquakes. One of the intriguing questions is whether similar 
scaling laws exist for other nonequilibrium systems. Another important puzzle is the origin of the unified scaling 
law. In the present paper we address the origin of pdfs in driven nonequilibrium systems, and we derive universal 
scaling laws connecting the magnitudes, recurrence-time, and spatial intervals of large events. These scaling laws 
should be applicable to all the above systems and other driven nonequilibrium systems. Here we also derive universal 
relationships between the scaling exponents. 

In the above mentioned driven-nonequilibrium system the large events do not have a typical energy scale. Hence 
under steady-state, the distribution of magnitude of large events follows a poweriaw. That is, 

P(I) ~ I~ a , 

where I is the strength of the events, and a > 0. Note however that the events cannot have energy more than what 
is available in the system, hence the distribution of the tail is decaying, mostly exponential. 

For temporal correlation study, we take events larger than a threshold value I c . Driven nonequilibrium systems 
can be classified in two classes, ones in which the events are uncorrelated, and the others in which the events are 
correlated. In the uncorrelated case, the pdf for recurrence of large events is Poissonian because the probability of 
large event is small. If the large events occur with certain rate R, then 1/R provides the time-scale for the recurrence 
of large events; consequently the pdf will be P(r) ~ i?exp(— Rt). The above feature is observed in BTW model 
[4|. In the correlated case, the nonequilibrium system is "critical" under steady-state, and the system does not have 
any time-scale for the recurrence time. Hence, the pdf of recurrence time is a poweriaw. We can combine the pdf of 
magnitude and recurrence time using the following law: 

P(/ C) r)=/-V- /3 /(-r/r c ) ! (4) 

where f(x) is an universal function. For small x, f(x) = 1, and for large x, f(x) — > exp(— x). Further application of 
scaling yields r c ~ IJ. It is expected that when .Li > I C 2, P(Ici, t) < P{I C 2, t), but r c i > t c i |7J,|ljj]. This observation 
implies that 7 > 0. Corral's law for earthquake @,la] is a good example which illustrates the above argument. 
The exponents a, (3, and 7 are not all independent. Using conservation of probability 

J P(I c ,r)dT = l, (5) 

and equating the powers of I c on both sides, we obtain 

a = 7 (l-/3). (6) 

As an example, the earthquake exponents of Corral |tJ, |ljj] are a rs 2/3, 7 « 1, /3 w 0.3, which satisfy the above 
relationship to a good approximation. 

Unfortunately Eq. © does not converge for P(I C , t) with (3 > 1 [Eq. <@}]. However we can circumvent this difficulty 
by computing the higher-order moments of r. If the integer part of f3 is n, then we compute 

{r n )= J r n P{I Cl r)dr. (7) 



Clearly the above integral converges. By postulating (t™) ~ J* n , and equating the exponents of I c on both sides, we 
obtain 

a = 1 {n+l-P)-5 n . (8) 

Hence the scaling arguments provide us an expression for P(I c ,t) [Eq. <@J] and relationships between exponents. 

The form of pdf of spatial separations between consecutive large events can be derived in the similar lines. The 
events could be either spatially uncorrelated or correlated. For uncorrelated case, the pdf of spatial separation is 
expected to be Poissonian: P(Ar) ~ exp(— Ar/(Ar) c )/(Ar) c , where the spatial density of the large events provide a 
length-scale (Ar) c . Nonequilibrium systems with spatially correlated events are critical under steady-state, and they 
do not have any length-scale for the separation between the consecutive large events. For this class of systems, the 
pdf is 

P(I C , Ar) = I- a (Ar)- c exp(-Ar/(Ar) c ), (9) 

with (Ar) c ~ Ic 1 *. If we assume dynamical scaling r ~ (Ar) z , then P(I C , Ar) = P(I C , r)dr / 'd(Ar) yields 

C = z(/3-l) + l. (10) 

Currently we do not know all the above parameters for any driven nonequilibrium systems. It will be interesting to 
test this conjecture. 

It should be noted that earlier Maslov et al. (jjji Paczuski et al. [ll|, and others had derived relationships between 
scaling exponents for several driven nonequilibrium systems, e. g., Bak-Sneppen model, Directed Percolation etc. 
[ill ]. These relationships differ significantly from our equations, which are based on large events. 

The scaling laws discussed above are quite general, and are expected to hold for systems with correlated large 
events. Earthquake is one such system for which Corral @,L!J| , Bak et al. [9J], and Davidsen and Paczuski [l(j have 
computed some of the scaling exponents. Solar flares, Bak-Sneppen model for punctuated evolution, and turbulence 
are other nonequilibrium systems which too have similar behaviour. It should be noted however that the exponents 
vary for different systems, but the form of pdfs is the same for all them. In the following discussion we apply the 
above scaling ansatz to ion-channel currents in voltage-dependent anion channels (VDAC), and test if it works. 

Our experimental setup on ion-channel current consists of two chambers containing buffer solutions KCl, MgCl<z, 
and HEPES. The chambers are separated by a thin wall of a polystyrene cuvette. A matched pair of Ag/AgCl 
electrodes, connected to Axopatch-200 amplifier, are put into the buffer solution of two chambers. A bilayer lipid 
membrane is formed on the pore. The protein molecules of VDAC purified from rat brain mitochondria (De Pinto 
et al. [20j) are put into the buffer solution of the outer chamber, after which they are stirred magnetically. The 
VDAC molecules now reconstitute the lipid membrane, and a channel is formed (Roos et al. |2l]]). After one channel 
is formed, the buffer solution (containing more protein molecules) of the outer chamber is replaced by fresh buffer 
solution. Thus we form a single ion-channel. 

Now the voltage is applied across the membrane, and the resulting current is recorded using the data acquisition 
software Clampex. This method is same as that used by Bezrukov et al. [22] and Banerjee et al [12] . In our experiment 
we measure currents with sampling frequency 1000 Hz for various applied voltages in the range — 30mV to 30mV. 
Typically the current has open and closed states with a fractal structure J23|. We find the length of a single open 
state to be too small. Therefore, we patch four long open-states and use the resulting timeseries for our analysis (see 
Fig. QJ • Due to a large variation of currents in various patches we perform the running average of the timeseries over 
eleven neighboring points. The noise is obtained by filtering the running average from the time-series. The statistics 
is performed on the resulting noise time-series. 

We take the absolute value of the current fluctuations, and compute the pdf of the magnitude of current fluctuation 
P(I). The log-log plot of P(I) vs. I is shown in Fig. The slope of the plot provides us the exponent a = 3.86±0.08. 

For the computation of temporal correlation, we filter the current using a threshold value I c . From the filtered 
current, we compute time- interval between subsequent current peaks. From this set of time-interval we compute the 
pdf of recurrence time P(I C , t). We perform this exercise for three values of I c : 3a, 3.5ct, 4ct, where a is the standard 
deviation of the current fluctuations. The corresponding plots P(I C , t) vs. r are shown in Fig. All the three plots 
have approximately the same exponent, (3 = 1.47 ±0.05. Note that for higher I c , P(I c ,r) is lower, but the powerlaw 
extends up to higher r. This feature is also seen in earthquakes. 

To obtain the universal function, we need to scale both x and y axes. If we scale only the y-axis to P(I, t)I"t^ and 
plot of P(I, t)I"tP vs. r, then we obtain a reasonable compactification along y axis, but not along x axis. See Fig. 



Q]for illustration. When we scale x-axis by r c = 12 and plot P(I,t)I"t@vs. rl~ a (see Fig. EJ, the compactification 
is quite good. The resulting plot is the universal function f{x) of Eq. (£jj. Note that we have taken 7 = a, which 
is an approximation. The actual parameters 7 and <5i [Eq. ||Hl] are not accessible in the present experiment. Since 
the current flow in the ion-channels is a localized, we do not have the pdf for spatial separation between consecutive 
events. 

In summary, the pdf of magnitudes of events in a driven nonequilibrium systems follow a powerlaw due to lack 
of any typical energy scale for the events. The pdf of recurrence-time could be either an exponential or a powerlaw 
depending on whether the events are temporally uncorrelated or correlated. The pdf of recurrence time for the 
many correlated nonequilibrium systems is a powerlaw due to lack of any timescale in the system. When the pdfs of 
magnitudes and recurrence-time for the large events are powerlaws, they could be combined into a single universal 
scaling law. In the present paper we derive the universal scaling law as well as relationships between the scaling 
exponents for driven-nonequilibrium systems. We show that universal scaling law is satisfied by ion-channel currents. 
Earlier, Corral P,Q and Bak et al. [9] had derived similar law for earthquakes. We have also derived similar scaling 
law for the pdf of spatial separation between subsequent large events. We should however keep in mind that the 
exponents vary for different systems, but the form of pdfs is the same for all them. 

Temporal and spatial correlations between events are seen in solar flares, turbulence, interface depinning models, 
music, punctuated equilibria etc. A timeseries of major historical events would presumably have temporal and spatial 
correlations. We should check how well scaling laws describe the above correlations. Higher-order correlations between 
the events should be another source for understanding nonequilibrium systems. 
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Figure Caption 

Figure 1: Timeseries of ion-channel currents in open state of Voltage-dependent Anion Channels (VDAC). Four 
parts have been patched here for extending the data range. 

Figure 2: Plot of the probability distribution of magnitudes of current fluctuation in ion-channel. The best fit is a 
powerlaw with slope a = 3.86 ± 0.08. 

Figure 3: Plot of the probability distribution of recurrence time P(I C , t) for three different threshold currents 
I c = 3a, 3.5a, 4a, where a is the standard deviations. The best fit for all the three plots is a powerlaw with the 
exponent (3 = 1.47 ± 0.05. 
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Figure 4: Plot of P(7 c ,t)t /3 vs. r. 

Figure 5: Plot of P(I c ,t)t^I" vs. Tl~ a . The resulting plot is the universal function f(x). 
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Figure 1: Timeseries of ion-channel currents in open state of Voltage-dependent Anion Channels (VDAC). Four parts have 
been patched here for extending the data range. 
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Figure 2: Plot of the probability distribution of magnitudes of current fluctuation P(I) in ion-channel. The best fit is a 
powerlaw with slope a — 3.86 ± 0.08. 




Figure 3: Plot of the probability distribution of recurrence time P(I C , r)for three different threshold currents I c = 3a, 3.5cr, 4a, 
where a is the standard deviations. The best fit for all the three plots is a powerlaw with the exponent (3 = 1.47 ± 0.05. 
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Figure 5. 



Figure 4: Plot of P(l c ,r)T (i vs. r. 
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Figure 6. 



Figure 5: Plot of P(I c ,t)t p I" vs. r/ c a . The resulting plot is the universal function f(x) 



